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Abstract. We present sub-second, continuous-coverage photometry of three flares 
on the dM3.5e star, EQ Peg A, using custom continuum filters with WHT/ULTRACAM. 
These data provide a new view of flare continuum emission, with each flare exhibiting 
a very distinct light curve morphology. The spectral shape of flare emission for the two 
large-amplitude flares is compared with synthetic ULTRACAM measurements taken 
from the spectra during the large 'megafiare' event on a similar type flare star. The 
white light shape during the impulsive phase of the EQ Peg flares is consistent with the 
range of colors derived from the megafiare continuum, which is known to contain a Hy- 
drogen recombination component and compact, blackbody-like components. Tentative 
evidence in the ULTRACAM photometry is found for an anti-correlation between the 
emission of these components. 



1. Introduction 

Active M dwarfs are notorious for white light flares in the near-UV and optical. White 
light emissio n is one of the first t ypes of radiation t o appear during both solar and 
stellar flares dNeidig & Kandll993l : lHawley et al.ll 1995b with a quickly-evolving impul- 
sive p hase followed by an extended decay phase that can sometimes last up to several 



hours (lHawley & Pettersenlll99ll : iKowalski et al.ll2010ah . Time-resolved spectral ob- 



servations of M dwarf flares show a rising continuum into t he near-UV, which has 
been modelled as a hot , T ~ 10,000 K blackbody spectrum (Ha wley & Fished 1 19921 : 



Fuhrmeister et all 12008). Recent spectral observations have revealed that a Hydrogen 
recombination component (B aimer continuu m) can also contribute significantly to the 
decay phase of the white light for A < 3646A ( Kowalski et alj|2010ah . 



Spectral observations can be severely limited by time resolution, with integration 
times (minutes) typically much longer than the timescales of flare evolution (seconds); 
as a result, broadband filters are more often used to characterize the time-evolution of 
the white light continuum. An analysis of 4 flares on the dM3e star AD Leo using 
\JY+UBVR filters has shown a nearly isothermal T ~ 8500 - 10, 000 K component in 
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the white light at all times (lHawley et al.ll2003l ). Another study by lZhilyaev et all (l2007r) 



analyzed high time-cadence (~0.1s) UBVRI photometry and concluded that multiple 
components (a Hydrogen recombination component and ~ 18,000 K blackbody) con- 
tribute to varying degrees during a flare on the dM3.5e star EV Lac. 

The origin of the white light continuu m during flares is a long -outstanding prob- 
lem in solar and stellar flare physics. The lAllred et aU (l2006h non-LTE radiative hy- 
drodynamic (RHD) models of M dwarf flares predict two continuum sources, but the 
photosphere does not receive enough energy to generate the temperatures near ~ 10,000 
K implied by the observations. However, the model white light emission resembles the 
shape of a hot (T ~ 8900 K) black body when convolved with the UV+UBVR filters 
(see Figure 12 of lAllred et al1l2006l) . indicating that there might be a degeneracy in the 
measured flare shape when using broad filters that encompass numerous spectral lines 
and continua. 

We have begun a flare monitoring campaign with ULTRACAM using custom 
continuum filters to characterize the white light continuum properties on sub-second 
timescales. In these proceedings, we present our first ULTRACAM observations of the 
flare star EQ Peg A. 



2. Observations 

ULTRACAM is an ultrafast, triple-beam, cascade dichroic CCD camera, and has helped 
open up simultaneous multiband, sub-second time domain astronomy dDhillon et al. 



20071) . On UT 10 August 2008, we obtained flare observations of EQ Peg A and B 
using ULTRACAM on the 4.2-m William Herschel Telescope at La Palma. EQ Peg 
A and B are a resolved binary system consisting of dM3.5e and dM4.5e c omponents, 
which have well-studied flare properties from the X-ray through the opt ical (lLacy et all 
119761 : iRobrade et al]l200l iMathioudakis et al.ll2006l : iLiefke et al.ll2008l) . 

The observations were taken with an Hor filter (FWHM ~ 50A) in addition to two 
custom continuum filters, NBF3510 (A c ~ 3510A, FWHM ~ 100A) and NBF4170 
(A c ~ 4170A, FWHM ~ 50A), designed specifically to avoid line emissi on and the 



B aimer jump at A = 3646 A. Figure 1 (left panel) shows a flare spectrum from Ko walski et al 
(2010a) with the NBF3510, NBF4170, and Johnson U band transmission efficiencies. 
The U band contains high-order Hydrogen emission and straddles the Balmer jump 
whereas the custom filters effectively isolate the continuum emission. Exposure times 
were 0.291 s for the NBF3510 camera and 0.134 s for the NBF4170 and Ha cameras. 

The data were reduced with the ULTRACAM pipeline software^, and we used a 
nearby comparison star, Gl 896 C, to obtain differential photometric measurements in 
all bands. A fixed aperture radius of 6 pixels was used. The night was photometric with 
variable seeing, and the observations were terminated due to morning twilight. 



3. Light Curve Data 

During the last ~30 minutes of observations, EQ Peg A produced three flares, each 
with a unique light curve morphology (Figure 1, right panel): a 'gradual' flare (G) at 



'http://deneb.astro.warwick.ac.uk/phsaap/software/ultracam/html 
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t ~ 50 min, an 'impulsive' flare (I) at t ~ 64 min, and a 'traditional' flare (T) at t ~ 70 
min. The high-cadence, continuous coverage photometry allows us to analyze the time- 
evolution of these flares without the ambiguity imposed by integration and readout 
times, which can be much greater than the timescales of flare evolution. We speculate 
that the differences in these light curve morphologies could be due to different active 
region structures (i.e., small vs. large arcades of flaring loops), different non-thermal 
particle beam flux evolution, and/or different white light continuum properties. 

In general, each flare has similar morphologies in the NBF3510 and NBF4170 fil- 
ters, and Table 1 summarizes the key properties. The G flare contains several emission 
peaks and a decay time that is long for its small peak amplitude. The T flare has a typ- 
ical fast-rise, exponential-decay ('FRED') shape, but exhibits several variations from 
this trend. The I flare is nearly symmetric about the peak, but a low-amplitude, short 
decay is apparent when examined over a small time window. It is interesting that the T 
and I flares have similar peak amplitudes, but the T flare has nearly 5x the equivalent 
duration (energy) as the I flare. 




Wavelength (A) Time ( m i nutes ) - MJD 54688.1971703 

Figure 1. Left - E xample flare (bl a ck) and quiescent (purple; scaled by a fac- 
tor of 3) spectra from iKowalski et all ( 2010a) with the transmission curves of the 
NBF3510, NBF4170, and the lMaiz ApeMnizI d2006h Johnson U band. The Balmer 
jump wavelength is indicated by a vertical dotted line. Note that the flare-to- 
quiescent contrast is much greater in the NBF3510 filter. Right - Three flares on 
EQ Peg A in the NBF3510 and NBF4170 filters. We show eveiy 4th point (top) and 
every 8th point (bottom) for clarity. 



4. Flare Color Analysis 

If the underlying spectral energy distribution of the star is known, then the relative count 
light curves of Figure 1 (right panel) can be converted to a light curve of the spectral 
shape of the flare emission. Flux-calibrated M dwarf spectra near the atmospheric cut- 
off are rare, so we obtained deep 3400 - 5500A spectra of EQ Peg A & B with the ARC 
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Table 1 . NBF35 1 Flare Properties 



Flare type 


peak amplitude (<x) 


Decay Time Constant [min] 


Equivalent Duration [min] 


gradual (G) 


2.1 (0.1) 


14 


7.6 


impulsive (I) 


5.1 (0.3) 


1 


3.2 


traditional (T) 


4.3 (0.3) 


7 


14.8 



3.5-m at the Apache Point Observatory in October, 2008. The spectra were reduced 
and flux-calibrat ed using stand ard IRAF0 procedures. Followin g the spectrophotome- 
try procedures of ISirianni et al. 1 (120051) and lMafz Apellanizl (12006). we estimate the ratio 



of quiescent fluxes of EQ Peg A in the NBF35 10 and NBF4170 filters to be 0.40, with 
~5% uncertainty. If C(t) is the relative count flux (normalized to 1 during quiescence), 
C is the relative count flux immediately before the flare, and Rq is the quiescent spec- 
tral ratio (^3510,2/^4170,2), then the spectral shape of flare emission around the Balmer 
jump is 

n C35io(?) ~~ Co ,3510 

RqX : . (1) 

Q170W _ C 

o,4170 



This formula follows from the integrand of the equivalent duration (IGershbergl 119721) 
for two filters, and we refer to it as the. flare color. 

4.1. The Flare Color as a Proxy for Two Continuum Components 

Recent spectral observations of a 'megaflare' on the dM4.5e star, YZ CMi, revealed an 
anti-correlated time-evolution betwee n the Balmer continuum (BaC) and a T ~ 10, 000 
K blackbody continuum components ( Kowalski et alj|2010al ). which c an be explained 



by a B almer jump in absorption forming during the secondary flares dKowalski et all 
l2010bl) . Figure 2 shows how this flare would appear in the ULTRACAM contin- 
uum filters, using synthetic photometry and flare colors calculated from the spectra. 
The NBF3510 emission contains both blackbody and BaC components, whereas the 
NBF4170 contains primaril y blackbody emis sion (and also probably a small amount 
of Paschen continuum; see lAllred et all (l2006h ). The anti-correlation between the total 
flare emission (U band, synthetic NBF3510, synthetic NBF4170) and the flare color is 
prominent throughout the entire decay and reflects the changing contributions of the 
blackbody and BaC components. 

For fast time cadences not achievable with spectral measurements, the ULTRA- 
CAM flare color acts as a proxy for the relative contributions of the two white light 
continuum components. The flare color evolution is shown for the I and T flares in 
Figure 3. The errors are obtained by propagating the pipeline uncertainties, the stan- 
dard deviation of the relative count flux in the preflare level, and the uncertainty in the 
quiescent spectral shape. The low-level flux enhancements during most of the flaring 
times lead to large errors for the original cadence. Binning the flare color into A? = 3.8 
s bins allows for a time-resolved analysis during the flare peak times. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation 
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Figure 2. The U band light curve for 1.3 hrs of the megafiare on the dM4.5e star 
YZ CMi (see lKowalski et al] feOlOal) for the complete light curve). From the simul- 
taneous optical ARC 3.5-m spectra, we calculated synthetic NBF3510 and NBF4170 
light curves and the flare color. An anti-correlation is prominent between the relative 
count fluxes and the color evolution. Note also the overall decline in flare color. 



The flare colors of the I and T flares have several notable similarities to the 
megafiare. The average flare color around the main impulsive phase of the I and T 
flares is ~2.2, which falls within the range of flare colors derived from the megafiare 
spectra, ~1.5 - 2.5. The I and T flares have much smaller amplitudes than the megafiare, 
which was emitting at very high levels, ~40 times the quiescent level, during the spec- 
tral observations. However, the similar flare colors imply a common geometrical scal- 
ing of the continuum components near the peaks of the I and T flares and during the 
decay phase of the megafiare, with the blackbody source having ~1/10 the area of the 
BaC-emitting region. 

We find tentative evidence for an anti-correlation between the flare color and the 
NBF3510 light curve at several times during the I and T flares. This is especially ev- 
ident during the fast rise and fast decline phases of the I flare, which exhibit similar 
anti-correlated variations as in the secondary flares at t ~ 95 min and ~ 130 min in 
Figure 2. If all continuum components varied by the same factor (i.e., if the increasing 
flare flux resulted from scaled up emitting area of all components), then the flare color 
would be flat as a function of time. The changing flare color likely indicates the vary- 
ing contributions of each continuum component, such that the blackbody has its largest 
relative contribution during the peak phase. Whether or not the decrease in flare color 
also results from the apparent amount of BaC emission decreasing as in the megafiare - 
due to increased B aimer continuum in absorption from the blackbody-like component 
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Figure 3. A closer view of the NBF35 10 data for the T flare (top panel) and I flare 
(bottom panel) reveals a variety of substructure beyond a simple fast-rise, exponen- 
tial decay shape. At this time-resolution, the rise in both flares has several distinct 
phases. The flare color is shown in red (right axis) and is binned by 3.8 s to reduce 
the errors. Note that the T flare has a decay phase that extends past the end of the 
observations at t ~ 83 min. 
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( Kowalski et alj|201 0b) - is difficult to determine from two-band photometry. Nonethe- 



less, the observed flare color of ~2.2 is consistent with this scenario. 

Despite the similarities between the I and T flares and the megaflare, one can see 
that the flare color in Figure 2 exhibits an overall decrease, whereas a decreasing flare 
color is not found in the decay of the I and T flares. This inconsistency is not yet fully 
understood, but it may be related to the presence of several secondary (blackbody-like) 
flares during the megaflare decay. 

4.2. Comparison to Radiative Hydrodynamic Flare Models 



We ca lculated NBF3510 and NBF4170 fluxes of the Fl 1 flare spectrum from the lAllred et al 



(2006) IHD models and found that the model flare fluxes and colors are inconsistent 



with the observed and derived properties of the I and T flares. Given an extremely large 
flare that emits from 5% of the visible hemisphere (i.e., more than twice the inferred 
area of the BaC-emitting region during the megaflare), the model predicts relative peak 
fluxes that are only 2.0 and 1.07 in the NBF3510 and NBF4170 filters, respectively. 
The discrepancy is particularly noticeable for the amplitudes of the I and T flares. The 
G flare has an amplitude that is more consistent with the model; however, the required 
surface area coverage is likely not realistic for such a short-lived, weak flare. The flare 
color of the Fll model spectrum is ~5, which is also inconsistent with the observations. 
Increasing the area of the ba ckwarmed photo sphere, which might be a more realistic 
geometry for a flaring region (llsobe et al.ll2007 ), would help to decrease the flare color 



of the model from such a discrepant value. 

It is a known problem that the hot blackbody component of the white light contin- 
uum is not produced in the RHD flare models. These models may fall short in matching 
the observed properties of flares due to not having enough power in accelerated elec- 
trons, and/or they may be limited by the one-dimensional treatment of a process that 
often includes arcades of flare loops (see th e spatial evolution of solar flare white light 
footpoints in Wan g et alJl2007l : IWandl2009h . 



5. Conclusions 

We use ULTRACAM continuum observations of EQ Peg A to characterize white light 
flares with custom continuum filters that do not contain Hydrogen line contamination. 
A variety of flare light curve morphologies are observed on sub-second timescales, with 
all of these flares exhibiting substructure beyond the standard fast-rise, exponential 
decay (FRED) shape. The rate at which each flare shape occurs is not known, but 
such information would help to understand the physical parameters and conditions that 
cause, for example, an I-type flare to occur instead of a T-type flare. 

We calculate a flare continuum color to quantify the spectral shape around the 
Balmer jump. The flare color evolution shows preliminary evidence of variations on 
timescales of several seconds during the impulsive phase. The average color over the 
first peak of these two EQ Peg flares is consistent with the synthetic flare colors during 
the decay phase of a megaflare on the M dwarf YZ CMi, for which we have spectra and 
therefore a detailed understanding of each continuum component. We find tentative 
evidence for an anti-correlated time-evolution between the relative flux and the flare 
color, similar to that previously presented for the YZ CMi flare. 

These ULTRACAM observations show that it is possible to use the NBF3510 and 
NBF4170 filters as proxies to understand the components of the white light continuum 
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on timescales not accessible by spectral observations. Additional observations with 
ULTRACAM of larger flares (hence, well-measured flare colors) would provide a more 
precise determination of the flare color during all phases of the flare. More observa- 
tions are also critical for finding and characterizing the properties of the anti-correlation 
between continuum components. 
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